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We investigate the efTects of impurities, either correlated disorder or a single Gaussian defect, on 
the collective dipole motion of a Bose-Einstein condensate of ^Li in an optical trap. We find that 
this motion is damped at a rate dependent on the impurity strength, condensate center-of-mass 
velocity, and interatomic interactions. Damping in the Thomas-Fermi regime depends universally 
on the disordered potential strength scaled to the condensate chemical potential and the condensate 
velocity scaled to the speed of sound. The damping rate is comparatively small in the weakly 
interacting regime, and in this case, is accompanied by strong condensate fragmentation. In situ 
and time-of-flight images of the atomic cloud provide evidence that this fragmentation is driven by 
dark soliton formation. 
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I. INTRODUCTION 

The creation of Bose-Einstein condensates (BECs) of 
ultracold atomic gases [H-Q has enabled investigations 
of some of the most fundamental concepts of condensed 
matter physics One of the most fruitful avenues of 
research has involved the use of BECs to probe the nature 
of superfluidity itself. Early studies led to observations 
of the critical velocity for the onset of dissipation [B-0| 
and quantized vortices [8l-[ll| . 

Recently, there has been much interest in using BECs 
to emulate disordered superfluids (c.f. Results 
from such experiments have wide ranging implications, 
from the transport of superfluid He in porous media 
to the motion of atomic BECs in microchip traps or mat- 
ter waveguides [isl - fisj . Of particular interest is how dis- 
order can disrupt, or even destroy, superfluidity. Due 
to their exquisite controllability, atomic BECs arc ideal 
physical systems with which to systematically study the 
interplay between superfluidity, disorder, and interatomic 
interactions. 

In this paper, we report measurements of the dissipa- 
tion of the superfluid flow of an elongated BEC subject to 
cither a disordered potential or a single Gaussian defect. 
We characterize the superfluid nature of the harmonically 
trapped cloud through detailed measurements of the ve- 
locity dependent damping of the collective dipole mode. 
We use a BEC of ""Li in the \F = l,mF ~ 1) internal 
state, where the interactions may be tuned via a wide 
Feshbach resonance located at 737 G [l9l - l2]| . This res- 
onance includes a shallow zero-crossing that enables the 
s-wave scattering length a to be tuned over a range of 
nearly 7 decades, with a as small as 0.01 oq, where oq is 
the Bohr radius [21[ . The gas may be made nearly ideal 
with transport properties strikingly different from the 
more strongly interacting case. Furthermore, the heal- 
ing length ^ = 1/ ^/Snnoa, where no is the peak density 
of the condensate, may be made as large as the conden- 
sate itself. In this regime, effects due to the fundamental 
wave nature of individual atoms become important. For 
example, if ^ is on the order of the disorder grain size 



or larger, a BEC can become an Anderson localized in- 
sulator [23, 1131 • In addition, the chemical potential fj, in 
this weakly interacting regime may be less than the ra- 
dial harmonic oscillator ground state energy, effectively 
"freezing out" the radial dynamics and leading to quasi- 
onc-dimensional behavior. 



A. Superfluidity of a BEC 

One of the seminal results originating from the theory 
of superfluid ^He is Landau's criterion. According to 
this criterion, elementary excitations can be created only 
if the fluid velocity v is greater than Landau's critical 
velocity wl [Hii] 



vl = nim . 

P 



(1) 



where e{p) is the energy of an elementary excitation of 
momentum p. For the case of a weakly interacting BEC 
with uniform density n, Bogoliubov theory gives the ex- 
citation energy as [26| 



(2) 




where m is the atomic mass and c is the bulk speed of 
sound. For small p, this spectrum reduces to the well 
known relation e(p) = cp describing phonon excitations 
with 



(3) 



where U = Airfi^a/m. Application of Eq.[T] gives vl = c, 
implying that only supersonic flow can dissipate en- 
ergy through the creation of elementary excitations; con- 
versely, if the flow is subsonic, excitations are energeti- 
cally forbidden, and the flow is superfluid. Application of 
Eq. [1] to the case of a non- interacting condensate implies 
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that vl = 0, suggesting that superfluidity cannot exist 
in an ideal gas. 

The dynamics of highly elongated BECs can be accu- 
rately modeled using an effective one-dimensional (ID) 
nonlinear Schrodinger equation (NLSE) [2^^ - [30| . In such 
a treatment, one starts from the 3D Gross-Pitaevskii 
equation (GPE), and integrates out the radial dimension. 
The effect of this integration is a reduction in c relative 
to Eq. [3] due to the average over the nonuniform radial 
density. For the case of a harmonically trapped BEG in 
the Thomas-Fermi regime, the bulk density n is replaced 
with the average density no/2. Therefore, the speed of 
sound becomes 

A theoretical description of an elongated BEG beyond 
the standard ID NLSE leads to a reduction (on the or- 
der of 10%) to the speed of sound relative to Eq. |4] [sij . 
In addition, the spectrum of axially propagating excita- 
tions in a cylindrical BEG can differ dramatically from 
Eq. [5] when /i ^ fiLOr, where ^ is the chemical potential, 
leading to an additional reduction vhvl (Hj. The highest 
/i condensates created in our system have fi/huj^ ^ 13, 
resulting in a predicted 20% reduction [3^ . 

When attempting to explain the onset of dissipation 
in any particular experimental situation, care must be 
taken to apply Landau's criterion locally^ by using the 
local density n{r ~ 0, z) instead of uq = n{r = 0, z = 0) 
in Eq. |4] |33| . For arbitrary trapping potentials, exci- 
tations will be nucleated first in regions of low density 
where the local speed of sound is small, and the critical 
velocity is reduced relative to the bulk. As a consequence 
of this effect, experimentally observed critical velocities 
arc often much lower than the bulk speed of sound [11-01 ■ 

The remainder of this paper is organized as follows: In 
Sec. In] we describe our experimental methods for creating 
a BEG in either a disordered harmonic potential or a har- 
monic potential with a single Gaussian defect; in Secs. lIIll 
and IIVI we discuss our results for the induced dissipa- 
tion for these two scenarios, where both the 3D Thomas- 
Fcrmi and the quasi-lD weakly interacting regimes are 
discussed for each case. We conclude in Sec.|V]with a dis- 
cussion relating the similarities and differences between 
dissipation in the two types of potentials, and directions 
for future studies. 



II. EXPERIMENTAL METHOD 

We create a BEG of ""Li in a highly elongated, cylindri- 
cally symmetric, hybrid magnetic-optical dipole trap pll . 
[s^l with radial and axial trapping frequencies in the 
ranges of a;^/ (27r) - 220-460Hz and w^/ (27r) ~ 4- 
5.5 Hz, respectively. The radial confinement is dominated 
by the optical trapping potential formed by a single fo- 
cused laser beam with wavelength 1030 nm and a 1/e^ 
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FIG. 1. (Color online) Disordered potential created from laser 
speckle, (a) Cut through an image of the speckle potential. 
The disorder strength Vd is proportional to the average value 
of the intensity (I) (dashed line), (b) The autocorrelation 
of the intensity distribution is well fit by a Gaussian with 
l/e^ radius (jd ~ 5.5 ^m. For some of the data in this paper 
(Figs. nil and [TOl fro = SA/im. 



Gaussian radius of 33 /im, while the axial confinement is 
dominated by an adjustable, harmonically confining mag- 
netic field. A set of Helmholtz coils provides a uniform 
bias field along the long (z) axis of the trap, allowin g fo r 
the tuning of a via a Feshbach resonance at 737 G [13- 
lU. The BEG is created at a field of 717 G where a is 
positive and large enough (~200 qq) to allow for efficient 
evaporative cooling in the optical trap, but small enough 
to avoid substantial three-body losses. At this field, the 
trap lifetime is limited to ~10 s due to three-body recom- 
bination with a loss coefficient of ^3 ^ 10~^^ cm^/s [35| . 
After evaporation, the BEG has no discernible thermal 
component from which we estimate that the temperature 
T < 0.5 Tc, where Tc is the BEG transition temperature. 
The bias field is then ramped over a timescale on the or- 
der of seconds to achieve the desired value for a. 

We excite the collective dipole mode of the conden- 
sate by pulsing on an axially oriented magnetic gradient, 
thereby abruptly shifting the center of the harmonic trap. 
After 1/4 of an oscillation period, the condensate is at the 
peak of an oscillation and we abruptly switch on either a 
disordered potential with an extent exceeding the oscil- 
lation amplitude of the condensate, or a single, narrow 
Gaussian defect located near the trap center. By vary- 
ing the duration of the gradient pulse we precisely vary 
the amplitude A of the oscillation, and therefore the ini- 
tial peak velocity vq of the condensate center of mass, 
where vq = Aoj^- At various times thereafter we image 
the cloud to track the center of mass location as well 
as the shape of the density distribution. We investigate 
the dependence of the damped dipole motion on vq, the 
strength of the disordered potential or single Gaussian 
defect, and on the value of a. 

The disordered potential is an optical speckle pattern 
created by passing a laser beam through a diffuser plate 
in a manner similar to previous studies [13, [36l - [38j . This 
beam is directed perpendicular to the trap z-axis. Fig- 
ure [T] shows a characteristic intensity slice of the disorder. 
The disorder speckle size ctd is defined to be the 1/e^ 
radius of a Gaussian fit to the autocorrelation of the in- 
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tensity pattern and is measured to be au = 5.5 fj,m. The 
beam has been cyhndrically focused such that in the ra- 
dial direction the speckle size is much larger than the 
radial Thomas- Fermi radius ^10 /im, making the disor- 
der effectively ID. We have verified that the intensity 
distribution of the disorder follows a decaying exponen- 
tial P{I) = {r}~^e~^/^^\ as expected for fully devel- 
oped speckle [S^. The average value of the speckle in- 
tensity (/) determines the disorder strength through the 
relation Vb = ?i (/) /(4/sat A), where the transition 
linewidth F = {2tt) 5.9 MHz and the saturation intensity 
/sat = 5.1 mW/cm^ The detuning from the '^Li 25" 2P 
transition A = (27r) 300 GHz, producing a repulsive dis- 
order potential. For the strongest disorder used in these 
studies, off-resonant scattering from the disorder occurs 
at a rate of ^ 0.1 s~^. The statistical properties of the 
speckle pattern are measured by direct imaging with a 
CCD camera before the optical system is installed onto 
the experimental apparatus. 

A cylindrically focused laser beam is used for the 
studies involving a single Gaussian defect. This 
beam has a Gaussian intensity distribution I{z,r) — 
loe^^'^^ /«'r+2 /fj)^ with beam waists Wr = 5mm and 
Wz = 12 /im. The radial size of the defect Wr is much 
larger than i?TF, ensuring that flow around the defect 
is suppressed. We conduct experiments using both a re- 
pulsive (blue detuned) and an attractive (red detuned) 
defect with |A| = 300 GHz. 

We adjust the healing length through an approximate 
range 0.5 ^j.m < ^ < 20 ^j.m by tuning a. Thus, a wide 
range of values arc achievable for the relevant dimension- 
less quantities, 0.1 < ^/ctd < 3.6 and 0.04 < < 1.7. 

III. DISORDER INDUCED DISSIPATION 

A. Thomas-Fermi Regime 

Figure [2] shows the position of the center of a con- 
densate at various times during a dipole oscillation in 
a disordered potential. The dipole oscillation is initi- 
ated by a kick that produces an initial peak velocity 
of vq = 20mm/s when the condensate passes through 
the center of the trap. For this data, the condensate 
begins its motion well into the supersonic regime with 
Vq 4 Co . The resulting oscillation is characterized by 
a time-dependent damping, suggesting that the damping 
depends on vq- The damping rate is initially small, goes 
through a maximum after about 3.8 s, and then dimin- 
ishes at later times. We fit 4-period sections of the data 
in Fig. [3] to the form of a damped harmonic oscillator: 

z(t) = ^e-'^*cos(w't + 0), (5) 

where tu' = (lo^ — The peak velocity vq is then 

computed from the fitted A for each data subset to ob- 
tain the damping coefficient /? as a function of Vq, with 
the results shown in Fig. [3] The damping monotonically 
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FIG. 2. (Color online) Damping of a condensate initially 
traveling supersonically through a disordered potential with 
Vo/h = 280 Hz. The center of the BEC (circles) is extracted 
from a Thomas-Fermi fit to the radially integrated column 
density (the "axial density"). The thick lines tracing the am- 
plitude are phenomenological guides to the eye. The initial 
amplitude is yl = 0.6 mm yielding an initial peak velocity 
of vo = 20mm/s. For this data, ujz = (27r)5.5Hz, ijjr = 
(27r)260Hz, a = 25 ao, and fj. = ^mujlR^^ = ft (1.1 kHz), 
where _Rtf is the axial Thomas-Fermi radius. In addition, 
Co = 5.6mm/s, ^ = 0.8 /im, and ^/ctd = 0.2. The insets show 
details of the oscillation at early and late times. 
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FIG. 3. (Color online) Velocity dependent damping. Results 
of fitting the data of Fig.[2]to Eq. [5]using a traveling 4-period 
window. The peak velocity no is obtained from uo = Aljz. 
The solid line is a square-root function convolved with an 
exponential decay and is meant as a guide to the eye. The 
inset shows the same data on a semi-log plot, emphasizing 
the nearly exponential decay of /3/a;- for large vo/cq. Vertical 
error bars correspond to the range in /3 for which A^^ = 1 for 
the fit to Eq. [5] while simultaneously adjusting A, /3, and (p to 
minimize x^- Horizontal error bars are determined using an 
identical process for A in Eq. [5]and are typically ~15%. The 
effects of systematic uncertainty in imaging magnification and 
variations in N are ~10% in the horizontal axis and ^5% in 
the vertical axis, these are not included in the displayed error 
bars. 
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increases for small vq , peaking near vq ^ 1 . 1 cq , followed 
by a nearly exponential decay of /3 for vq > Cq. 

A perturbative theoretical treatment has produced a 
closed form solution for the velocity-dependent damp- 
ing, resulting in good quantitative agreement with our 
measurements [40|. For weak disorder the qualitative 
behavior shown in Fig. |3] can be understood through a 
local Landau critical velocity argument. At low veloci- 
ties, Bogoliubov quasiparticles are only created within a 
thin shell at the surface of the condensate, where the low 
density leads to a low local speed of sound, and there- 
fore a low local vl- As the velocity of the condensate 
increases, a larger condensate volume can support exci- 
tations because a larger fraction of the atoms violate the 
local Landau criterion. The maximum damping occurs 
near the point where the velocity of the BEC reaches the 
peak speed of sound cq in the condensate. At even larger 
velocities the excitation volume cannot increase further, 
but the Bogoliubov density of states decreases, resulting 
in the observed exponential decrease of the damping. 

Except for the absence of a critical velocity, the quali- 
tative behavior of the velocity dependent damping shown 
in Fig. 13] is remarkably similar to that predicted by ID 
NLSE simulations of a uniform, repulsive BEC in the 
presence of an oscillating Gaussian obstacle [4l|, |42| . In 
these simulations, above a certain impurity strength- 
dependent critical velocity, the impurity moving at a ve- 
locity V deposits energy into the BEC in the form of den- 
sity fluctuations. The average rate of condensate energy 
growth (dE/dt) increases nearly linearly with v, to a peak 
at ?j ^ c as the defect excites dark solitons and linear 
sound waves. As the velocity of the defect is increased 
further, the density fluctuations decrease signiflcantly, 
accompanied by an exponential decrease of (dE/dt), sim- 
ilar to our experimental observations. In contrast to a 
single impurity in a uniform condensate, a defect is al- 
ways present in a low density region of a condensate in 
a disordered harmonic trap. Consequently, vq is always 
greater than the local speed of sound at the edge of the 
condensate and excitations are always present. Previous 
experimental js^jHEIill and numerical [11] studies of the 
damping of collective modes and the damping of Bloch 
oscillations in a disordered lattice potential [45,|46| have 
found qualitatively similar results. 

Figure S] shows in situ polarization phase-contrast im- 
ages [13] of the BEC at various times in the oscillation 
shown in Fig. [2] The damping clearly does not result 
from a loss of collectivity as predicted by ID NLSE nu- 
merical simulations H^]. Rather, the BEC nearly main- 
tains its original shape throughout the oscillation. Close 
inspection of the density distributions in Fig. |4] reveals 
a "tail" of non-condensed atoms that appears to oscil- 
late slightly out-of-phase with the central Thomas-Fermi 
distribution. At early times, these non-condensed atoms 
appear to lag behind the BEC, while at later times they 
oscillate in-phase with it. This two-component out-of- 
phase oscillation is reminiscent of the second sound-like 
oscillation reported in Ref. . In that work, the initial 




FIG. 4. (Color online) Characteristic in situ polarization 
phase-contrast images of the data shown in Fig. [5] at various 
times. The images are nearly equally spaced in time between 
the time labels. 
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FIG. 5. (Color online) Generation of a non-condensed com- 
ponent, (a) Squares show the center of the Thomas- Fermi 
(condensed) component and circles show the center of the 
Gaussian (non-condensed) component. The Gaussian cen- 
ter trails behind the Thomas-Fermi center and has a smaller 
amplitude of oscillation. Within experimental uncertainty, 
ujz ~ (27r) 5.1(2) Hz for both components. For this data, 
a = 200ao, = 3 x 10^ n/h = 1.8kHz, Vu/f- = 0.22, 
vo = 28mm/s, co = 7.2mm/s, and Ur = (27r)220Hz. (b-d) 
Axial density distributions with bimodal fits (solid lines) and 
a single component Thomas-Fermi fit (dashed lines) at var- 
ious times during the oscillation: (b) 28 ms, (c) 100 ms, (d) 
190 ms. The condensates in (b) and (d) are traveling in the 
positive direction whereas the condensate in (c) is traveling 
in the negative direction. 
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Disorder strength V^/^ 

FIG. 6. (Color online) Damping vs. Vd- Open circles corre- 
spond to the data shown in Fig. [7] (a = 200 ao) in the range 
0.7 < vo/co < 0.9; filled circles correspond to a = 200ao, 
vo/co = 1.2, fj,/h = 2.2kHz; squares correspond to a = 25 ao, 
vo/co ~ 1.2, ^i/h = 750 Hz. The damping parameter /3 follows 
a power law with p ~ 2 (solid and dashed lines) , independent 
of ^ or a. To minimize systematic effects associated with the 
velocity dependence of /3 (e.g., Figs.[2]and[3|, we fit a 4-period 
window for which the data is described well by Eq. (5] Verti- 
cal error bars are as defined in Fig. O The inset shows the fit 
values of p as a function of vq / co for a collection of data sets 
at a = 200 ao. The dashed lines indicate the plus-and-minus 
one standard deviation extent for the collection of measured 
velocities. Vertical error bars for p are determined as in Fig. |3] 
using a fit to Eq. [S]for each oscillation at a given vq/cq. Data 
corresponding to filled circles and squares was taken using 
an optical trap setup different from that described in Sec. [H] 
with A = 1064 nm and a beam waist of 24fj,m resulting in 
to^ = (27r)4.9Hz, LOr = (27r)460Hz, and iV = 3 x 10^ Also, 
for these data sets ao = 3.4 ^m. 



temperature was high enough that damping occurred due 
to the interaction between a BEC and a thermal compo- 
nent. In contrast to those results, we observe that the 
dipole oscillation is undamped in the absence of the dis- 
ordered potential. Furthermore, there is no observable 
heating due to the quick switch on of the disorder. In our 
experiment, therefore, the presence of the non-condensed 
component seems to be linked to the motion of the BEC 
in the disordered potential. A recent numerical simula- 
tion using a truncated Wigner method predicts the emis- 
sion of incoherent atoms from a BEC moving supcrsoni- 
cally through a disordered potential |49[ , consistent with 
our observations. 

We have investigated this effect in further detail using 
in situ absorption imaging, which allows for determina- 
tion of the low density non-condensed wings of the dis- 
tributions. Figure [5] shows that by fitting the cloud to 
a bimodal Thomas-Fermi plus Gaussian profile, a phase 
difference of At/) = 0.23 between the condensed and non- 
condensed cloud centers is found. Note that the interac- 
tion strength is different for this data than that shown in 
Figs. 

We have systematically investigated the dependence 
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FIG. 7. (Color) Transport regimes of a BEC traveling through 
a disordered harmonic potential. Black squares show the val- 
ues of disorder strength Vo/i-i and initial peak center of mass 
velocity vq/cq for the data used to extract /? from a fit to Eq.[5] 
using 4-6 periods of oscillation. The interpolated color map 
(and contour lines) for P/cOz is derived from the measured 
results. These measurements have a — 200 ao, N — 2 x 10^ 
atoms, fi/h = 1.5 kHz, cq = 6.5mm/s,ajr ~ (27r)260Hz, and 
ujz = (27r) 5.5 Hz. The variable experimental quantities are A 
and Vd. Due to small shot-to-shot fluctuations in the posi- 
tion of the center of mass of the cloud, measurements with 
vo < 0.2 Co are not reliable. Data with /3 < 2 x 10"'^ is con- 
sistent with undamped motion. 



of /3 on the disorder strength Vd . Figure [S] shows the 
normalized damping parameter I3/uJz plotted against the 
normalized disorder strength Vd/Mj where /i is the chem- 
ical potential of the condensate prior to the kick and 
before the disorder is switched on. We find the data fits 
well to a power law 



/3 



oc 



(6) 



for all measured velocities. The precise value of p, how- 
ever, depends weakly on vq across the range of velocities 
< vq/co < 5, with a mean value oi p = 2.1(5) (see 
Fig. [H] inset). 

Figure [7] presents the measured values of /3 as a func- 
tion of both Vd and vo- As expected, a vertical trace 
through this plot shows a qualitative similarity to Fig. |3l 
We observe two distinct regimes of reduced damping: one 
where vq/cq <^ 1 and the other when vq/cq 3> 1, with 
the damping reaching a maximum at wq ~ cq. A numeri- 
cal simulation using an effective ID NLSE has produced 
qualitatively similar results [i^ . 



B. Variation with Interaction Strength 

Wc observe nearly universal behavior for /3 as a func- 
tion of both Vd//!, and vq/cq for BECs in the Thomas- 
Fermi regime. As already shown in Fig. [6l /3 oc (Vd //i)^ 
for condensates with fi differing by a factor of 3. Shown 
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FIG. 8. (Color online) Universal damping vs. vq/cq. The 
disorder strength was adjusted to keep 0.30 < Vo/n < 0.35 
for all of the data. Squares correspond to a — 28 ao, — 
2.5 X 10^ fj,/h = 550 Hz, co = 4.0mm/s, uj^ = (27r)5.5Hz, 
and Ldr = (27r) 260 Hz; open circles correspond to a = 200 ao, 
TV = 3x 10^ fi/h = 2.4kHz, cq = 8.3mm/s, = (27r)4.5Hz, 
and iUr — (27r)460Hz; filled circles correspond to the same 
parameters as Fig. [T] Error bars are as defined in Fig. [31 
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FIG. 9. (Color online) Peak damping vs. a with fixed Vd/a* 
and vo/cu. For this data, Vb and vq were adjusted to keep 
0.3 < Vd/m < 0.4 and 0.6 < vo/co < 1.4 with ah other pa- 
rameters as in Fig. [T] The upper horizontal axis shows values 
for fi obtained from a variational solution of the GPE [2l|] 
(note that the upper tick marks are not strictly logarithmi- 
cally spaced). The linear fit has a slope 0.002 O;^^. Vertical 
error bars are as defined in Fig. El 



in Fig. [5| is a comparison between the damping at inter- 
action strengths a = 200 ao and a ~ 28 ao, with con- 
stant Vu/fJ.- Although the respective values of cq differ 
by nearly a factor of 2 between the two data sets, the 
peak damping occurs at ?;o/co ~ 1 for both, demonstrat- 
ing the nearly universal behavior of /3 vs. vq/cq. On the 
other hand, the peak damping rate between the two data 
sets differs by nearly a factor of 5, showing that while 
the general shape of the damping curve is universal, the 
magnitude of the damping is not. 

An investigation of the effect of interatomic interac- 
tions on the peak damping {vq/cq ^ 1) at fixed Vj^/fi 
is shown in Fig. |9| We find that /3 scales linearly with 
a, going to zero with decreasing interactions, consistent 
with the disappearance of the low energy phonon portion 
of the excitation spectrum as U — > 0. 

The elongated confinement geometry in our system fa- 
cilitates the investigation of the dimensional crossover 
from the 3D to the quasi-lD regime where fi <^ TiuJr 
[sol [sH ]. Shown in Fig. (TOjare measurements of /3 vs. /i 
at constant and vq. When fi > huir (to the right 
of the vertical dashed line) and vq is comparable to, or 
less than cq (as is the case for the data shown as open 
and closed circles), we find /3 oc /i^^-^. By reference 
to Fig. [71 one can gain a qualitative understanding of 
this behavior going from high to low fi: starting subson- 
ically (open and filled circles) , the system travels along a 
path from the weakly damped regime (lower left corner 
of Fig. [7]) towards the regime of strong damping (middle 
right region). As fj, decreases, the quantities Vn/fi and 
Vq I Cq increase correspondingly, and the system follows a 
path which crosses several contours of constant /3 while 
approaching the strongly damped regime near w/c ~ 1. 
Consequently, the system displays a strong dependence 
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FIG. 10. (Color online) Damping vs. /i with fixed Vd and 
Vq. Squares and filled circles correspond to Vd//i = 370 Hz 
with Wo = llmm/s and 6mm/s, respectively. Open circles 
correspond to Vo/ft = f40Hz and vq = 6mm/s. The vertical 
dashed line denotes \i = Tiujr, at which point dq/co = 1.7 for 
the open and filled circles and vo/co = 3 for the squares. We 
varied ^ by adjusting a, shown on the upper horizontal axis 
(note that the upper tick marks are not strictly logarithmi- 
cally spaced). Values for ^ are obtained from a variational 
solution of the GPE [2lt| using the following measured exper- 
imental parameters: oJr = (27r)460Hz, ujz — (27r)4.5Hz, and 
= 4 X 10^ atoms. For this data an = 3.4 /im. Vertical error 
bars are as defined in Fig. [3l 



of l3 on ^. Blue squares depict a different situation where 
the system is supersonic for all /i investigated. For large 
fi the system occupies a point in Fig. [7] with vq/cq > 1 
and Vd/m < 1- As decreases, the system follows a di- 
agonal path, roughly tracing a contour line of constant 
/3, moving into the regime of vq/cq 3> 1 and T^d/m 3> 1 
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FIG. 11. (Color online) Damping of a nearly non-interacting 
gas. For this data a = 0.4 ao, iV = 2 x 10^, fi/h = 26 Hz, c = 
1.2mm/s, Vo = 4: fj., ujr ^ (27r)240Hz, and oj;, = (27r)5.3Hz. 
(a) Center of mass position (circles) and radius (squares) of 
the condensate as a function of time. Here we use statistically 
determined values for the center of mass Zcm = J z n{z)dz / N 
and radius R, given by i?2 = 4/(2 - z^mf n{z)dz/N. (b-d) 
Axial density traces at various times in the oscillation: (b) 
40 ms (c) 100 ms, (d) 960 ms. After two full oscillations, the 
cloud has fragmented and spread to a size comparable with 
the initial oscillation amplitude. 



the collective dipole mode brought on by extensive frag- 
mentation. Because > /x, it is perhaps not surprising 
that the condensate quickly fragments. While the cen- 
ter of mass of the cloud damps after about 5 oscillation 
periods, examination of shot-to-shot differences in the 
damped density distributions reveal that the position of 
the fragments are highly non-rcpcatable, suggesting that 
some fragments remain in motion. This residual motion 
is consistent with the long thermalization time expected 
from weak two-body interactions. It is interesting to 
note that the maximum single particle kinetic energy, 
Ek = ^mijjlA^ = h (295 Hz), is 2.8 times larger than the 
average height of the disordered potential. The observed 
dephasing is therefore consistent with the expected be- 
havior of a gas of non-interacting particles interacting 
with a disordered potential where the disorder strength 
is smaller than the kinetic energy of the individual par- 
ticles. 



IV. DISSIPATION INDUCED BY A SINGLE 
GAUSSIAN DEFECT 

A. Thomas-Fermi Regime 

In an effort to better understand the mechanisms re- 
sponsible for the damping by disorder, we have investi- 
gated the dissipation induced by a single Gaussian defect. 
The defect potential is described by V{z) — Vue'^^ Z™^, 
where Wz = 12 /xm. The static effect of either an at- 
tractive or repulsive defect on a repulsively interacting 
BEG in the Thomas-Fermi regime is shown in Fig. [T^l 
As expected, the attractive defect leads to an increase of 
the density in the region of the defect, accompanied by 
a small decrease of the density in the wings of the distri- 
bution, while the opposite is true for a repulsive defect. 



(top-right corner of Fig. [7]). When /i < hui,. (to the left 
of the vertical dashed line), and uo 3> cq, we observe a 
negligible dependence of /? on fi. In this quasi-lD regime, 
/3 is affected only by changing Vu or wq: consistent with 
the behavior expected for a nearly ideal, classical fluid. 
This may be understood by reference to Eq. [2] where for 
V ^ c, the first term in the Bogoliubov excitation spec- 
trum dominates making the system "quasi-ideal" with 
e{p) independent of fi. 

Figure [TT] shows damping of a weakly interacting gas 
with a = 0.4 ao, deep into the quasi-lD regime, where 
n/hutr ^ 0.1. We find that Vd produces the 

same damping {P/uJz = 0.07) as that for a BEG with 
a = 200 ao and Vu = 0.25 fi. The nature of the damped 
motion of a weakly interacting gas in strong disorder 
is strikingly different from the damped motion of a 
strongly interacting gas in weak disorder, even though 
the timescale of the damping in both cases is compara- 
ble. Figure [TT] shows that the damping in the weakly 
interacting regime is caused by the loss of coherence of 



The dynamical distributions can differ dramatically 
from the static shown in Fig. [T3] where in situ 

axial densities are displayed for various times through- 
out the dipole oscillation. In the following discussion 
we refer to the upstream side of the condensate as the 
portion of the BEG that reaches the barrier after the 
leading or downstream portion. The interaction of the 
BEG with the repulsive defect, shown in Fig. [T^b). pro- 
duces a deep downstream density rarefication as well as a 
large upstream density compression bearing a qualitative 
similarity to a shock wave. Similar structures have been 
predicted in effective ID theoretical treatments [13, [13 
and interpreted as upstream and downstream dispersive 
shocks. In contrast, the interaction of the BEG with the 
attractive defect, shown in Fig. [T^ a). produces no such 
shock waves. However, the cloud is slightly compressed 
near the defect simply due to the attractive defect po- 
tential. Because vq > cq, phonon excitations cannot be 
emitted in the upstream direction as they would have to 
propagate faster than the speed of sound. Glose inspec- 
tion of Fig. I13f a) reveals minimal density modulation of 
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FIG. 12. (Color online) BEC in a harmonic trap with a single 
Gaussian defect, (a-c) correspond to a repulsive defect, while 
(d-f) to an attractive one. (a) and (d) in situ polarization 
phase-contrast images, (b) and (e) axial densities correspond- 
ing to the images, and (c) and (f) numerical solutions to the 
GPE with the dashed lines showing the solution in absence 
of a defect. The inset trace shows the characteristic shape 
of the potential. For all panels, a = 200 ao, A*' = 4 x 10^, 
uj:: = (27r) 5.0 Hz, and ujr = {2tt) 360 Hz. 
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FIG. 13. (Color online) Axial densities at various times dur- 
ing a supersonic oscillation in the presence of an (a) attrac- 
tive or (b) repulsive defect. The arrows are proportional to 
the instantaneous velocity of the condensate. The vertical 
dashed lines denote location of the defect. For this data, 
a — 200 ao and = 13 mm /s; (a) corresponds to A'^ = 4x10^, 
/i = 1.5 kHz, vo/ca = 2, Vd = -0.8 ^i, cj^ = (27r)4.7Hz, 
and ujr = (27r)360Hz; (b) corresponds to A" = 1 x 10*^, 
fi = 3kHz, vo/co = 1.4, Vd = 0.4^, uj^ = (27r)5.0Hz, and 
ujr = (27r)360Hz. 



the upstream side, while more modulation is evident on 
the downstream side. 

Several ID theoretical studies predict the formation of 
downstream propagating dark solitons in addition to an 
upstream dispersive shock as a repulsive defect is super- 
sonically swept through a condensate (sol . 13, [53456| , 
which is consistent with the density fluctuations visible 
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FIG. 14. (Color online) Density fluctuations produced cross- 
ing a repulsive defect. Absorption images after 4 ms time-of- 
flight. The vertical dashed line denotes the location of the de- 
fect. The experimental parameters are as stated in Fig. llSr b'). 



in Fig. [131 However, the size of the dark solitons will 
be on the order of the healing length, ^ = 0.5 fim for 
these condensates, which is a factor of 6 smaller than our 
imaging resolution. 

Figure [M] shows time-of-flight images of the BEC os- 
cillating in the presence of a single repulsive defect. In 
contrast to the in situ images of Fig.fTHlb'). after time-of- 
flight additional structures emerge which were not previ- 
ously visible. These structures are consistent with dark 
solitons that form from short length scale in situ phase 
fluctuations that map onto larger scale density modula- 
tions after time-of-flight. However, in situ phase fluctu- 
ations may also arise from thermal excitations in highly 
elongated BECs, and these can also manifest as density 
fluctuations after time-of-flight [13 • Close inspection of 
Fig.HHreveals that deep density modulations are present 
only in the downstream portion of the BEC (after the 
flrst pass through the defect), consistent with the dark 
soliton interpretation. Similar density fluctuations have 
also been interpreted as dark solitons in an experiment 
using a moving defect and a stationary BEC [58j . 

We have measured /3 as a function of both Vd and 
vo, with characteristic results presented in Fig. [13 Con- 
trary to what was observed for a disordered potential, 
we observe a critical velocity Vc below which the motion 
is undamped, for both the attractive and repulsive de- 
fects. We flnd that the peak damping for an attractive 
defect is signiflcantly weaker than for a repulsive one. 
Figure [Tni presents measurement results of /3 as a func- 
tion of both Vd and vq. For an attractive defect, we flnd 
that Vc/cq ~ 0.7 with Vc depending only weakly on Vd. 
However, for a moderately strong repulsive defect, Vc/cq 
occurs signiflcantly below 1 and depends strongly on V^d • 
For both attractive and repulsive defects Vc tends to cq 
as IVd/a*! is reduced to zero. 

Once again, a model based on a local Landau criterion 
is sufflcient to explain the dependence of Vc on Vd- For 
simplicity, consider a uniform density flow impinging on 
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FIG. 15. (Color online) Velocity dependence of /3 induced by 
a single Gaussian defect, (a) Attractive defect with (squares) 
Vd/^J = —0.8 or (circles) Vo/n = —0.3 and other parame- 
ters as stated in Fig. llSf al except with N = 8 x 10^ and 
fj,/h = 2kHz. (b) Repulsive defect with (squares) Vd/p = 0.4 
or (circles) Vu/fJ- = 0.2 and other parameters as stated in 
Fig. Iisr b). Both types of impurities show critical behavior 
at low velocities as well as undamped motion at large vq/cq- 
Note the difference in scale between damping induced by an 
attractive versus a repulsive impurity. Vertical and horizontal 
error bars are as described in Fig. |3l 
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FIG. 16. (Color) Transport regimes of a BEG traveling 
through a harmonic potential with a central Gaussian de- 
fect. Coordinates of the black squares are the values ofVu/fJ, 
and wo/cq for the data sets used to extract /S/ojz from a fit to 
Eq.O The color map (and contour lines) for P/uj^ is derived 
from an interpolation using the measured results. Dashed 
white lines show the local Landau critical velocity as given by 
Eqs. [7]and[9l The attractive and repulsive cases are qualita- 
tively similar: superfluidity for vo/co ^ 1, increased damping 
as vo/co — >■ 1, and reduced damping for wo/cq 3> 1. Damping 
induced by an attractive impurity is an order of magnitude 
weaker than for a repulsive one. Data with Vd < and Vd > 
correspond to parameters in Figs. llSf a') and (b), respectively. 



either a repulsive or attractive Gaussian potential [5J]. 
With the assumption that the superfluid flow pattern is 
stationary, the local density of the condensate near the 
defect must be modified in a similar way to that shown 
in Fig. [T^for a static defect. For the repulsive case, the 
local density is reduced near the defect, resulting in a 
lower local speed of sound. In addition, flux conservation 



requires that the local condensate velocity increase in the 
low density region near the repulsive impurity to preserve 
the stationary flow pattern. A corresponding argument 
can be made for the case of an attractive defect. These 
effects serve to increase the local value of v{z)/c{z) near a 
repulsive defect and decrease it for an attractive one. As 
a result, excitations can be created near the peak of the 
repulsive defect in a BEG with a center of mass velocity 
that is significantly lower than the bulk speed of sound. 
For the case of an attractive impurity, on the other hand, 
one expects excitations to occur in the bulk condensate 
first, rather than near the impurity, and therefore at a 
flow velocity near the bulk speed of sound, as observed. 

We quantify this picture, in the case of a repulsive de- 
fect, by applying the local Landau criterion at the instant 
the center of the BEG crosses the defect. In Ref. 4^ the 
authors used an effective ID NLSE in the high density 
regime to determine the locus of points where the local 
condensate velocity v{z) is equal to the local speed of 
sound c(z); this defines the curve 



Co 



= 1- 



5/2 



yo > 0, 



(7) 



where Vb //.t = Sno /no is the fractional change in the 
peak density at the peak of the repulsive defect. When 
vq/cq < 1, we can ignore effects of the axial Thomas- 
Fermi profile of the condensate because A ^ Rtf for 
our trap, where Rtf is the axial Thomas-Fermi radius. 
Equation [7] is plotted in Fig. [12] when Vb > and is 
found to agree with the measured Vc for the range of Vb 
explored experimentally. Therefore, the observed reduc- 
tion of the critical velocity below cq is consistent with 
the local Landau critical velocity without invoking more 
exotic mechanisms, such as vortex nucleation. This is in 
contrast with several experiments involving BEGs in less 
elongated configurations @, 0: as well as in superfluid 
^He where nucleation of vortex lines and rings can result 
in Vc < vl [25| . 

In the case of the attractive defect, the density, and 
therefore c{z), is enhanced at the location of the defect 
and reduced only slightly elsewhere. We find that the re- 
duction in density in the bulk due to the enhancement at 
the defect is less than 1% for the strongest barriers used, 
leading to an essentially unperturbed speed of sound in 
the bulk. The ratio of the local fiuid velocity to the local 
speed of sound can then be found by considering only the 
bare Thomas-Fermi profile, and is given by 



c{z) Co \l~z^/R'^ 



^TF 



(8) 



where, using Eq. IH vo/co = 2 A/ Rtf- If 2 A < Rtf 
then Vol Co < 1 and the local Landau criterion is satisfied 
everywhere inside the condensate: 



- = l,Vb < 0, 
Co 



(9) 



implying that Vc is independent of lb- Our measure- 
ments, however, show that Vc depends weakly on Vb with 
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Vc/cq 1 only in the weak impurity limit. Our experi- 
mental results are consistent with numerical simulations 
using a ID NLSE [i^] for which the local Landau crite- 
rion accurately describes the repulsive impurity case, but 
slightly overestimates Vc in the attractive case. 

Figures [15] and [1^] demonstrate that damping is sig- 
nificantly suppressed deep into the supersonic regime. 
We observe undamped motion when vq is greater than 
a Vn-dependent upper critical velocity Numerical 
simulations [12, 144. l59l| have shown that for "wide and 
smooth" barriers «C Wz) the emission of radiation from 
the defect in the form of phonons and solitons can be 
very small for supersonic velocities. In fact, it has been 
shown analytically that the radiation emission rate re- 
sulting from a defect moving supersonically through a 
condensate decreases exponentially with the ratio £,/wz 
[60| . Without emission of radiation, energy dissipation 
is inhibited and the flow persists, even though Landau's 
criterion is violated. For the data presented in Figs. [12]- 
I17[ S,/wz ^ 0.04, well within the regime where supersonic 
non-dissipativc flow is predicted. Experiments similar to 
ours have also shown a reduction in soliton emission from 
a barrier moving through a condensate in the supersonic 
regime [sst . 

We therefore observe three distinct regimes of flow in 
the single defect system: subsonic superfluid {vq/cq < 1), 
dissipative (wq/cq ^ 1), and supersonic non-dissipative 
{vo/cq > Figure [17] displays axial densities from 

in situ polarization phase-contrast images at the instant 
the defect passes through the peak of the condensate for 
the three different velocity regimes. As expected, for the 
superfluid flow regime the axial density profiles look very 
much like the equilibrium profiles of Fig. 1121 there is 
an increase (decrease) in the density at the location of 
the attractive (repulsive) defect. In the dissipative flow 
regime, on the other hand, the fiow patterns for Vb > 
show significant distortion, while for Vd < there is 
little distortion, as discussed in detail above. Finally, in 
the supersonic non-dissipative flow regime, we observe 
a counter-intuitive density inversion with respect to the 
superfluid regime, where the attractive defect produces 
a density depression while the repulsive defect causes a 
density peak. 

The physical origin of this counter-intuitive density in- 
version can be understood by considering the behavior of 
the gas at large vq. In this regime, as in the disordered 
case, the Bogoliubov excitation spectrum, given by Eq.[2l 
is dominated by the /2m term and therefore, domi- 
nated by plane waves with wavenumber k = p/fi rather 
than phonons. For this "quasi-ideal" gas, the drag should 
be determined by the scattering of these plane waves off 
of the defect [SJ]. At high velocities, scattering of these 
waves from the defect is greatly suppressed, leading to 
low dissipation. If we extend this argument further and 
consider the atoms to be classical particles, one expects 
the atoms to slow down in the presence of the repulsive 
defect, resulting in a density increase near the defect, 
while the opposite is expected for an attractive defect. 
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FIG. 17. (Color online) In situ density distributions of a 
condensate passing through a Gaussian defect. These images 
are taken at the instant the center of the BEG first crosses 
the defect. Rows correspond to the three fiow regimes: sub- 
sonic superfiuid (vq/co < 1), dissipative {vo/co ~ 1), and 
supersonic non-dissipative {vo/co > 1). (a-c) Attractive de- 
fect with Vd/ij. ~ —0.85 and vq/cq = 0.31, vq/cq = 1.0, 
and vo/co = 3.0, respectively, (d-f) Repulsive defect with 
Vd/ fJ- = 0.65 and vq/cq = 0.15, vq/cq — 0.90, and vo/co = 2.0, 
respectively. The arrows indicate the direction and relative 
speed of the condensate. For this data, all other parameters 
are as described in Fig. [TS] 



Density inversions similar to the ones presented here 
have also been discussed in the context of dissipationless 
stationary states at supersonic velocities [H, [60l - l62j as 
well as sonic black holes [6^. Under our experimental 
conditions, when vq/cq ~ 1 the edge of the barrier can 
serve as a sonic event horizon. Such systems have been 
proposed as possible candidates with which to study "ta- 
ble top" astrophysics, where exotic effects, such as Hawk- 
ing radiation, should be observable. Interestingly, in this 
system the experimenter plays the role of the so-called 
super-observer, having access to the regions both outside 
and inside the event horizon [63l-[66l. 



B. Weakly Interacting Regime 

Figure [TH] shows results of measurements of a weakly 
interacting condensate (a = 0.6 cq) oscillating in the 
presence of a repulsive defect. Under these conditions the 
condensate is in the quasi-lD regime with ^/TiiOr = 0.15. 
We find that the axial density profile of the condensate 
becomes increasingly modulated during the damped os- 
cillation, consistent with theory [42,[55|. We compute the 
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FIG. 18. (Color online) Oscillation of a weakly interacting 
BEC in the presence of a repulsive defect, with a — 0.6 ao, 
AT = 2.5 X 10^ n/h = 44 Hz, Vd/h = 0.8, lj, = (27r)4.7Hz, 
iOr = (27r)300Hz, and vo/co = 1.6. (a) Center of mass po- 
sition as a function of time (computed as in Fig. Illf) giving 
P/oJz = 0.03; (b) Root-mean-square density deviation from a 
Gaussian fit to the axial density distribution (see text). The 
solid line shows the decay of the oscillation energy (in arbi- 
trary units) found from the fit in (a), (c-e) In situ axial den- 
sity traces and Gaussian fits at various oscillation times: (c) 
ms; (d) 140 ms, at the second crossing of the defect; and (e) 
1260 ms, after several crossings of the defect. At large times 
we find that the large density modulations are accompanied 
by only a slight increase of the axial size of the condensate. 



root-mcan-square deviation of the axial density distribu- 
tion from a Gaussian fit rifit as a proxy for the increased 
internal energy of the condensate due to the density mod- 
ulation 




n{z) - n^tiz) 



(10) 



where we take the integration length L to be over the 
central 70% of the condensate to minimize edge effects. 
Figure [TBl b) shows that A initially increases in time and 
then saturates. The time-dependent increase in the den- 
sity modulation qualitatively matches the loss of oscilla- 
tion energy. We therefore conclude that the damping of 
the dipole mode is caused primarily by the creation of in 
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FIG. 19. (Color online) Velocity dependence of damping 
with a repulsive defect. Circles correspond to a nearly 
non-interacting BEC with a = 0.6 ao, fi/h = 31Hz, co = 
0.9mm/s, and Vd ~ 0.9 /i. Shot-to-shot variations in the po- 
sition of the cloud limit the extraction of /3 to i;o > Imm/s, 
corresponding to do/co > 1.1. Squares correspond to data 
from Fig. I15f b) for comparison, a — 200 ao, At/ft = 3kHz, 
Co = 9.2mm/s, and Vd = OAfj,. Error bars are as defined in 
Fig. El 



situ density modulations in the cloud. 

Results of measurements of the velocity-dependence of 
the damping by a repulsive defect with a = 0.6 oq are 
shown in Fig. [121 As was the case with disordered po- 
tentials, we find that the timescale for damping in the 
quasi- ID regime with a strong impurity strength is much 
longer than that observed in the Thomas-Fermi regime 
with a weak impurity strength. 



C. Dark Soliton Production in the Weakly 
Interacting Regime 

Of particular interest in the quasi-lD regime is the 
ability to create and observe long-lived dark solitons. 
These nonlinear excitations have been previously created 
in BECs with repulsive interatomic interactions through 
a variety of means, including direct phase irnprinting 
[stI . [ssI. spatially selective microwave transfer j69| , slow 
light [70t , two condensate interference [7l|, [72| , and, sim- 
ilar to the work presented here, as a result of a BEC 
crossing a semi-permeable defect . 

In general, the decay of dark solitons occurs as a re- 
sult of dynamical instability or as a result of dissipative 
dynamics associated with the interaction of the soliton 
with quasiparticle excitations of the BEC. However, it 
is known that dark solitons can have very long lifetimes 
in the quasi-lD regime [73j . For the most weakly inter- 
acting BECs presented here, fx/hiOr = 0.13, making our 
system ideally suited to study long-lived dark solitons. 

We have studied the formation of deep in situ density 
modulations in BECs for different values of a, with the 
results shown in Fig. 1201 Dipole motion is initiated af- 
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FIG. 20. (Color online) Dark soliton formation. In situ ax- 
ial densities of BECs during the first (a-d) and fourth (e-h) 
passes through a semi-permeable defect. The defect is located 
at 2 = 0, and its strength was adjusted to keep Vd/a* ~ 0.7. 
Oscillation amplitudes were adjusted to keep vq ~ cq. (a), 
(e): a = 0.1 ao, A*' = 1.0 X 10^ fi/h = 5Hz, ^ = 12.5 ^m, 

= 16(6) Aim; (b), (f): a = 0.5 ao, = 2.2 x 10^ 
fi/h = 30 Hz, ^ = 4.9 Aim, = 6.6(2) Aim; (c), (g): a = 1.7ao, 
iV = 2.6xl0^ fi/h = 77 Hz, ^ = 3.06 Atm, = 2.8(4) Atm; (d), 
(h): a = 5.4ao, Af = 2.2 x 10^ fi/h = 144Hz, ^ = 2.24 Atm, 

— 2.5(3) Atm. The trap frequencies for this data are 
Ur = (27r)240Hz and cu^ = (27r)4.75Hz. The dashed lines 
show fits to Eq. 1111 We omit the fit in (d) for clarity. For 
comparison, the thin dashed line in (a) is only the Gaussian 
portion of the fit. Error bars for are given by the standard 
deviation of a collection of images. 



ter the field is slowly ramped to a desired value near the 
scattering length zero-crossing at 544 G. Panels (a-e) of 
Fig. [20] show the cloud after 3/4 of a complete oscilla- 
tion. The defect is switched off after the first pass of 
the cloud, and the cloud is imaged after it returns to the 
center of the trap after another quarter period. There- 
fore, ~ 100 ms elapses between the initial interaction of 
the cloud with the defect, where the soliton is created, 
and imaging. Deep density modulations, consistent with 
the formation of stable dark solitons, are observed. For 
comparison, panels (e-h) of Fig. show the cloud after 
passing through the defect 4 times. The density modula- 
tions in this case appear less monochromatic than in the 
single pass case, suggesting the presence of both linear 
(phonons) and nonlinear (solitons) excitations. We ex- 
tract the healing length by fitting the single-pass data 
in Fig.Hnito [zl 



n(z) = Ae 



1- D sech^ 



(11) 



where A is the background density, a is the size of the 
atomic cloud, D is the depth of the soliton, zq is the loca- 
tion of the soliton, and is the healing length. Through 
a variational solution of the GPE, we can independently 



estimate the healing length ^ using the measured values 
of N, a, LOz, and lOj.. The results of this analysis are re- 
ported in Fig. [501 The average size of the density dips is 
very nearly the healing length predicted by the GPE esti- 
mations, i.e. ^ ^. This observation is consistent with 
the formation of a downstream dispersive shock consist- 
ing of a train of dark solitons as a supersonic BEG crosses 
a semi-permeable barrier [H, IH, [53l - l56j . 



V. SUMMARY AND FUTURE DIRECTIONS 

We have conducted comprehensive measurements of 
the dissipation of superfluid flow in an elongated BEG 
subject to either a disordered potential or a single Gaus- 
sian defect. By measuring the velocity and disorder 
strength-dependent damping parameter, we have char- 
acterized the breakdown of superfluidity of a harmoni- 
cally trapped cloud in both the 3D Thomas-Fermi and 
the quasi- ID weakly interacting regimes. 

Our data largely support the validity of the Landau 
criterion for a critical velocity above which the super- 
fluid motion is damped, as long as the criterion is applied 
locally. The local criterion accounts for the inherent in- 
homogeneity of trapped gases, as well as density modifi- 
cations produced by large defects. The only exception is 
for attractive defects of relatively large strength, where 
we find that Vc decreases to Vc ^ 0.7 cq for Vu/ fJ, < —0.5. 
Dissipation is also found to diminish for velocities greater 
than which we associate with reduced excitation of 
dark solitons and phonons. 

Throughout the 3D Thomas-Fermi regime, the damp- 
ing is found to be well described by a universal relation 
depending on the dimensionless defect strength Vd / ^ 
and velocity vq/cq. The universal damping peaks at 
vq/cq ~ 1 for any Vd/a* and scales as (Vd/a^)^ for all fi. 
As At decreases, the peak damping rate decreases as well, 
consistent with the disappearance of the phonon portion 
of the excitation spectrum as cq — t- 0. Damping in the 
quasi-lD regime is qualitatively different. In this case, we 
find for fixed absolute Vd and vq that /? is independent 
of At. In this regime, damping is accompanied by frag- 
mentation and spreading of the cloud, with the damping 
monotonically increasing with Vd/Ek, where is the 
maximum single particle kinetic energy. 

A particularly intriguing possibility for the future is to 
explore the transport properties of a weakly attractive 
gas. In the case of a disordered potential there exists the 
opportunity to study the transport properties of bright 
matter-wave solitons with the prospect of observing 
Anderson localization in such systems [zl, [T^I- For a 
single defect, there is also a possibility for the creation 
of coherently split solitons or solitonic Schrodinger's cat 
states [zi-ig. 
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